ANALYSIS OF A STATE CHANGING SUPERSOFT X-RAY 

SOURCE IN M31 

B. Fate\\ R. Di Stefano^, T. Nelson^, F. A. Primini^, J. Liu^, S. Scoles^ 

ABSTRACT 

We report on observations of a SSS in M31, rl-25, that has exhibited spectral 
changes to harder X-ray states. We document these spectral changes. In addition, 
we show that they have important implications for modeling the source. Quasisoft 
states in a source that has been observed as an SSS represent a newly-discovered 
phenomenon. We show how such state changers could prove to be examples 
of unusual black hole or neutron star accretors. Future observations of this 
and other state changers can provide the information needed to determine the 
nature (s) of these intriguing new sources. 

Subject headings: galaxies: individual (M31) - X-rays: binaries 



1. INTRODUCTION 



Luminous supersoft X-ray sources (SSSs) were established as a class by Einstein Obser- 
vatory observa tions of roughly 30 sources in the Magellanic Clouds, Milky Way, and M31 
(IGreinerl 120001 ). Chandra and XMM- Newton observations of external galaxies have now dis- 
covered hundreds of soft X-ray sour ces with properties that both exemplify and extend the 
class of SSSs (e.g., lOrio et al.ll2010l . references therein). Even though they are bright, with 
luminosities higher than 10^^ ergs s~^, we know of only a handful in the Galaxy, because the 
radiation they emit is readily absorbed by the interstellar medium. In fact, the SSSs used to 
define the class display lit tle or no emiss ion above 1 keV. Roughly a dozen SSSs are known 
in the Magellanic Clouds (IGreinerl 120001). Some of these are associated with novae, and are 
clearly hot white dwarfs (see Greinei 12000 ). In M31, SSSs h ave been shown t o be associated 
with supernova remnants, novae, and active galactic nuclei (jOrio et al.ll2010l ). 
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The more mysterious component of the class is comprised of X-ray binaries, most 
with orbital periods of a day or less. A promising model for these sources is one in 
which their prodigi ous luminosities are produc e d by the nuclear burning of matter accreted 



by a white dwarf (Ivan den Heuvel et al.l Il992l : iRappaport et al.l Il994 



Hachisu et al. 1996 



Pi Stefano fc Nelson! Eggj h Nuclear burning should allow the white dwarf to retain ac- 
creted matter and increase in mass. Binary SSSs have therefor e been suggested as progen- 
i tors of accretion-induc ed collapse (Ivan den Heuvel et al.l 119921 ) and of Type la supernovae 
( iRappaport et al.lll994l ). Nuclear-burning white dwarfs in wider orbits are also expected . 
Indeed, symbiotic binaries have been observed as SSSs (e.g.. lGreinerll2000l : lOrio et al.ll2007l ). 



Because it is difficult to detect SSSs in the Milky Way, and the Magellanic Couds are 
too small to host a large population, it is important to search for SSSs in external galaxies. 
The advent of Chandra and XMM-Newton has made such searches possible. Hundred s of 
SSSs have now been discovered, some in galaxies as far from us as the Virgo cluster ( iLiu 
20081 ). As the numbers of SSSs has increased, we have begun to find evidence of sources that 
have properties different from those of the SSSs that established the class. Some SSSs are 
hundreds of times more luminous than the Eddington limit for a Chandrasekhar-mass white 
dwarf. Several of these ultra luminous supersoft sources are candidates for accreting black 
holes ( jPi Stefano et al.ll2004l ). In addition, the search for the softest sources has identified 
a class of sources that are significantly harder than SSSs, yet also significantly softer than 
canonical X-ray sources. Quasisoft X-ray sources h ave luminosities aljove 10^^ ergs s~^, but 



emit few or no photons with energy above 2 keV (jPi Stefano et al.ll2004[ ). Some could be 
hot white dwarfs in which there is an additional hard component and/or which are highly 
absorbed. Those fitting this model are good candidates for progenitors of Type la supernovae, 
because they would lik e ly correspond to the most massive nuclear-burning white dwarfs (e.g., 
Rappaport et al.lll994j : |Pi Stefanoll2010l . references). Others are too hot to be white dwarfs 
and may correspond to either black holes of neutron stars. 

In this paper we report on a source (rl- 2 5) that has been observed to switch bet ween 
SSS and QSS states Jstiele et al.lboosl . l20ld |Pi Stefano e^aPboiol : lOrio et al.l[201ol ). No 
such sources are known in the Galaxy or Magellanic Clouds. We examined all available 
Chandra and Swift data and optical data from HST and from the Local Group survey. We 
also checked the literature for XMM-Newton observations and analysis of the source. The 
question we want to answer is: what is the physical nature of this unique source? 

We analyze the X-ray and optical data for the source in §2 and §4, respectively. In §3, 
we present evidence that the source changes state. We discuss the possible models that fit 
this sources in 55. 
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X-ray Observations and Analysis 



The state changing source rl-25 has been reported on multiple times ( Kong et al.l 12002 



2010 



Di Stefano et al. 



2010 



WiUiams et al-lbooilPi Stefano et al.ll2 004l: l Kaaretll2002l IVoss fc GilfanovlboOTi : Istiele et al. 



Orio et al.. 2010i) . It is located in the central region of M31 



approximately 0.4' from the nucleus. The coordinates of the source are (J2000.0) RA = 
00:42:47.90, DEC = +41:15:49.99 (see §4 to see how it was calculated). This region of 
M31 has been well sampled over the past 11 years. In this paper, we searched the Chandra 
archive for all public observations of rl-25 through June 2, 2009. The observations that 
covered the source are shown in Table 1 and include 86 observations, rl-25 was detected in 
45 observations (28 ACIS-I, 2 ACIS-S, and 15 HRC-I observations) from August 8, 2000 to 
March 11, 2009. We note that the source position was the same in all detections, and we 
are confident that rl-25 is one source (see Figure 1). We supplemented the Chandra data 
with Swift observations of the region in 2009. However, the source was off during the Swift 
observations and are not included in this paper. 



We took photometry results from iLiul (120081 ) for observat ions l isted in that paper. We 



photometered the recent AC 



version 4.1.2 ( iFruscione et al. 



S dat a in the same manner as iLiul (120081 ). We used CIAO 
2006h to analyze the recent data. For s ource detection, we 



used CIAO tool wavdetect (IFreeman et al.ll2002t IFruscione et al.l 120061 ). We set the scales 
for wavdetect to 1", 2", 4", and 8". The chance detection threshold was set to 10~^. The 
source aperture size was defined as the 3(T elliptical source region computed by wavdetect. 
The background aperture was defined as an elliptical annulus around the source region. The 
annulus had radii of two and four times the source ellipse radius. Since the source is not near 
any other X-ray sources, the background aperture size was suitable for all ACIS observations. 

Only photons between 0.1-7 keV were considered and the following bands were used to 
classify them: Soft, 0.1-1.1 keV; Medium, 1.1-2 keV; and Hard, 2-7 keV. The counts in each 
band and the total for the source were computed as C^ — C;, x Ag/Ai,, and the uncertainty was 
computed as 1 + a/O.75 + + C;, x [As/Ai,Y, with Cg as the counts in the source aperture. 
As as the source aperture area, as the counts in the background aperture, and Ai^ as 
the background aperture area. The total count rate, as well as the rate in each band, were 
corrected by a vignetting factor. The vignetting factor was derived from the exposure map 
as the ratio between the local and the maximum map value. 



For HRC-I observations, we used the CIAO tool dmextract (IFruscione et al.ll2006l ) to 
extract raw counts in source and background apertures, as shown in Figure 2, and to com- 
pute net counts. We defined aperture sizes so that the source aperture encircled energy 
fr action was ~1 and t he background aperture fraction 0. We also used the CIAO aprates 



( IFruscione et al.l l2006l ) tool to compute the background- marginalized posterior probability 
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distribution for the source rates, assum ing non-informative prior distributions for source and 



background rates (IKashyap et al.ll2009l ). We note that the HRC data were not corrected for 



vignetting. Thus, the rates are lower than the ACIS counterparts by 10% or less. 

For sources with low flux significance, defined as netrate / anetrate < 3, we used the 
posterior probability distributions to estimate the 3a upper limits to the rate. The upper 
limits were calculated using aprates. The same method was applied to both HRC and ACIS 
observations to determine upper limits. For ObsID 10552, there were no counts in the source 
or background apertures. We could not calculate upper limits for this observation. 

Table 2 shows the photometry results for rl-25. The first column is the Modified Julian 
Date (MJD) for each observation. The next three columns are the counts in the Soft, 
Medium, and Hard bands, respectively. The fourth column is the Total counts for ACIS 
observations (in the energy range of 0.1-7 keV). The uncertainties are la in magnitude and 
calculated using Poisson statistics. The next four columns are the Soft Count Rate (SCR), 
Medium Count Rate (MCR), Hard Count Rate (HCR), and the Total Count Rate (TCR), 
respectively. The last column of Table 2 is discussed in §3.2. 

Figure 3 shows the light curve for rl-25. The top panel shows the count rate in each 
band (soft, medium, and hard) vs. obsdate, and the middle panel shows the total counts 
(in ACIS observations) vs. Obsdate. The top panel demonstrates that the relative count 
rates change over time for the source. The middle panel indicates that the source alternates 
between being dominated by soft and medium photons. The bottom panel of Figure 3 is 
discussed in 53.2 



3. Evidence for State Change 
3.1. XSPEC 

We used the CIAO tool dmextract to create source spectra for all observations where 
more than 60 source counts were collected. These were ObsIDs 1575 (ACIS-S, 183 source 
counts), 4720 (ACIS-I, 61 source counts), 4721 (ACIS-I, 65 source counts), and 4722 (ACIS-I, 
62 source counts). Source spectra were extracted from a circular region with radius 6 pixels 
centered on the source. Associated background spectra were extracted from an annular region 
(also centered on the source) with inner and outer radii of 10 and 25 pixels, respectively. 
For all observations, we crea. ted spectral response files using the CIAO tasks mkacisrmf and 



mkarf ( jFruscione et al.l 120061 ) . 



We modeled the resulting spectra using the XSPEC package (lArnaudlll996[ ). Since the 
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number of source coun ts is rather small, we chose to fit the unbinned spectra using the 
C-statistic ( lCashlll979l ) rather than the statistic, which has been shown to introduce a 
systematic bias in parameter estimation in low count rate spectra. One caveat in using the 
C-statistic is that the background is negligible within the source extraction region. This 
is definitely the case in our ACIS-I observations (~1 background count expected in the 
source region), and an acceptable approximation in the ACIS-S observation of rl-25 (only 
17 background counts). The C-statistic as implemented in XSPEC does not provide a 
measurement of the goodness of fit, but can be used to constrain parameter confidence 
intervals, which we have done here. 

We fit a simple absorbed blackbody model to each spec trum in the energy range 0.3-8. 
keV, using the wabs model for the interstellar absorption (IMorrison fc McCammonl Il983l ). 
None of the spectra had enough counts below 1 keV to place a tight constraint on the 
interstellar column density (in several cases the fits were consistent with no absorption), so 
we fixed A^^ at two values representative of the range expected towards a n X-ray source in 
M31 (1.1 X 10^^ and 6.4 x 10^^ cm~^). The lower limit was taken from Di Stefano et al. 



( I2OO4J ). The upper limit is shown for completeness; it is very unlikely that the source has 
such a large column density. The resulting parameter values and their 1 a uncertainties 
are presented in Table 3. We note that our results for ObsID 1575 (we found 0.120 keV 



< kT < 0.140) are consistent with those reported in iDi Stefano et al.l (|2004[ ) (they found 
that kT = 0.122 keV). Examining the temperatures found for rl-25, a clear increase in the 
later observations can be seen, independent of the choice of absorbing column. At the 90% 
confidence level, the temperature found for the ACIS-S spectrum is ~0.1 keV lower than 
the ACIS-I values. The blackbody temperature returned by the model fitting is constrained 
primarily by the high energy cut off. Although ACIS-I has poorer low energy sensitivity, the 
fact that significantly harder counts are detected in such short exposures indicates that the 
increase in model temperature is real and independent of the differences between the two 
ACIS instruments. 

In order to assess the statistical significance of our obser ved state chan ges, we used the 
simulation capabilities of the spectral fitting package XSPEC ( lArnaudlll996l ). In observation 
1575 (a 38 ks ACIS-S observation), 183 counts were detected for the source and we derived 
two best fit thermal models (depending on assumed Nh, see Table 4) with kT = 0.08 or 
0.13 keV. Later, in observations 4720, 4721 and 4722 (three 4 ks ACIS-I observations), 
approximately 60 counts were detected in each, and the derived temperatures were in the 
range 0.15-0.27 keV. We want to determine the likelihood that, given the differences in 
instrumental response and exposure time, the best fit model from ObsID 1575 produces the 
observed counts distribution in ObsIDs 4720-4722. 



- 6 - 



We used the fakeit command in XSPEC to simulate 1000 observations, using the 
best fit model from ObsID 1575 as input, and the appropriate response files and exposure 
time for the location of rl-25 in each of the ACIS-I observations. For each simulation, we 
recorded the counts in each of the energy bands outlined in §2, and then used the ensemble 
of simulations to determine the associated error. Finally, we compared this to the observed 
counts distribution. 

Our results are recorded in Table 4. It is clear that the softer best fit model found 
for observation 1575 cannot account for the observed counts distribution in observations 
4720-4722: such a source is simply too soft to give rise to the required counts above 1 keV. 
The same is true when considering the opposite situation — the best fit models derived for 
the harder state would produce a much larger number of counts in the ACIS-S band than 
observed, primarily because of the long duration of that particular observation. Looking at 
Table 4, the observed counts in all cases differ from the tested model by a significant margin 
(up to 20cr). Therefore, we can state with confidence that the source experienced a real 
transition from a softer state to a harder one. 



3.2. Spectral Estimates from pimms 



For the majority of ACIS detections of rl-25, there were fewer than 60 total counts. We 
found that using a spectral fitting prograr a was not a v iable option for these observations. 
Instead, we used HEASARC tool pimms (IMukail 119931 ) included in the CIAO package to 
estimate the appropriate spectral models. For cycles Chandra AOl and A02, we had to 
create effective area maps, as pimms no lo nger provides support for these cycles. The maps 



were created by running CIAO tool mkarf (IFruscione et al. 



20061 ) for on-axis sources in both 



cycles (note, mkarf was not run on rl-25, but for other X-ray sources that were on-axis). 
The process was completed for the ACIS-I, ACIS-S, and HRC-I instruments. Effective area 
maps created by mkarf differ slightly from the default pimms maps because they are created 
for specific observations. However, the different effective area maps produce similar spectral 
estimates, which we discuss later in the section. 



The method we used is similar to the iDi Stefano et al.l (120041 ) approach. The process 



MCR 



and 



HCR 



Tcw TCR^ ""-"^ TCR Same ratios calculated by pimms for 



compares the observed 

different spectral models (we call the pimms ratios y§^^ TCRp ^- pimms 

ratios are within the la uncertainty range for all three of the observed ratios, the model is a 
match. We used the same set of inputs for all ACIS observations that detected the source. 
Our XSPEC analysis suggested that black bodies fit the data well. So, we used a black body 
model to fit the data using pimms. The black body kT input was varied from 0.025 to 1.224 
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keV. We assumed an Nh — 1-1 x 10^^ cm~^. We input an unabsorbed X-ray flux (0.1-7 
keV) corresponding to Lx = 10^^ ergs s~^. This value was simply an initial guess for the L^; 
the guess was incorrect for the observations. However, we corrected the X-ray luminosity 
(calling it the Lx,p) using the ratio of the TCR and TCRp. That is: 

^x. = |^xl0^^er^..- (1) 

Table 5 shows an example of the pimms method for ACIS-S observation 1854 (OBSMJD 
51922.4). The first column is the kT (in keV) that was input. The next three columns show 
the ratios for the 'rcRp ^ TCRp ^ respectively. The last column shows the value of 

Lx for each value of kT, and were calculated using equation (1). The ratios in the table 
were compared to the observed ones (the values of /cT that fit the observation are in bold 
in Table 5). The range for the observed ratios were: 

SCR 

, MCR 
2.65 X lO"'^ < — — — < 0.262 
TO R 

MC R 

-8.43 X 10-^ < — — < 8.43 x 10"^ 
TO R 

We applied a second step to the method to reduce the range of temperatures that fit 
each observation. We broke the soft band down into smaller bins, as shown in Table 6, 
and followed the same procedure as above. That is, we compared the ratios of the observed 
values and those calculated by pimms for these smaller bins. This additional step significantly 
decreased the range of kT that fit each observation. For example, using just the first step 
for ObslD 1854 (sec Table 6), there were three values of A; T that fit (0.083, 0.112 kcV, 0.151 
keV). After applying the second step, we were able to exclude the two higher values of kT. 
However, for many observations the k T range was still over 100 eV after applying the second 
step. 

As a test of the pimms method, we applied it to the four observations that were analyzed 
with XSPEC. Table 7 shows a comparison of the fitting results from XSPEC and pimms for 
A^^ = 1.1 X lO^^cm"^. The first two columns are the observation ID and MJD. The third 
and fourth columns show the k T range determined by XSPEC and pimms, respectively. The 
two methods produce similar results, demonstrating the validity of our pimms estimation 
method. 
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Table 8 shows the results from the pimms method. The XSPEC results are shown in 
place of the pimms estimates for the four observations described in §3.1. The last two columns 
show the range of kT and Lx that fit each observation. For all observations, the minimum 
value of kT that fit corresponds to the maximum value of Lx- The table shows clear evidence 
for state change. This can be better seen in Figure 4, which graphically display the spectra 
for rl-25. Again, we used the XSPEC (black symbols) results instead of the those from 
pimms (blue symbols) for the four observations analyzed in §3.1. The plots show a change 
in both state and Lx for the source. The source is observed in a SSS state (at minimum k T 
< 100 eV), and changes to a QSS state (at maximum kT > 245). The source also exhibits 
a change in of ~1 magnitude. We note that, hke the XSPEC results, increasing the Nh 
decreases the kT and increases the Lx that fit the data. However, as with the XSPEC 
results, the change of state observed in the source using the pimms method is seen at larger 
values of N^. This implies that the state change is not dependent on the value of Nh chosen. 

As a check of the default effective area maps used by pimms, we created effective area 
maps for several observations. The effective area maps were created for on-axis sources 
using mkarf . We then recalculated the ratios predicted by pimms using the new maps. The 
ratios calculated from the two methods differed by 0.01 at most. This difference leads to a 
small change (at most 5 eV) in the range of kT values that fit. This is a negligible change 
considering the inherent uncertainties associated with using pimms. Furthermore, the kT 
range for most observations was ~100 eV (see Table 8). Thus, there would still be evidence 
for state change if mkarf effective area maps were used instead of the default pimms maps. 

We used the pimms estimates to convert ACIS-S and HRC-1 count rates to ACIS-I units 
(see the last column of Table 2, as well as the bottom panel in Figure 3). This calculation is 
done to compare the results between Chandra instruments. We choose ACIS-I units because 
most of the detections of the source were in this instrument. Thus, the last two columns of 
Table 2 show the same count rate for ACIS-I observations. 

For ACIS-S detections of rl-25, we use pimms to convert the total count rate to ACIS-I 
units using the appropriate spectral models (see Table 8). We ran pimms twice to determine 
the upper and lower bounds of the converted rate. To determine the upper bound, we input 
the following count rate: TCR + eTCR, where eTCR is the uncertainty in total count 
rate. The model used was the maximum value of kT determined to fit the observation. 
To determine the lower bound, we input: TCR — eTCR, using the minimum value of kT 
determined to fit the observation as the model. In both calculations, we assumed a Nh = 
1.1 X 10^1 cm-2. 

HRC-1 observations and ACIS-S non-detections of the source provided no spectral in- 
formation. For these observations, we assumed a spectral model and then applied the same 
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process as above to convert to an ACIS-I count rate. All observations in Chandra cycles 
AOl, A02, and AOS were converted assuming: 0.075 keV < kT < 0.200 keV. For the rest, we 
determined the appropriate range of kT from the first ACTS detections before and after the 
observation. For example, HRC-I observation 7285 (MJD 54052.3) detected the source. The 
first ACTS detection preceding it was ObsID 7140 (MJD 54002.8; 0.160 keV < fcT < 0.205 
keV, see Table 8) and the first after was ObsID 7064 (MJD 54073.9; 0.120 keV < A;r < 0.165 
keV, see Table 8). We used the following range for kT for ObsID 7285: 0.120 keV < kT 
< 0.205 keV. 



3.3. X-ray Color Ratios 



Using the method described in iPrestwich et al.l (120031 ) . we calculated the soft and hard 
color ratios for all rl-25 detections with 25 or more counts (the results are shown in Table 
9). The ratios were calculated using the following equations for each of the 12 observations 
in the table: 

{Hard Counts — Medium Counts) 



Hard Color Ratio 



Total Counts 



(2) 



Soft Color Ratio 



{Medium Counts — Soft Counts) 



(3) 



Total Counts 

The color ratio s were used to construct a X-ray color-color plot in the same fashion as 
Prestwich et al.l ( l2003l ) (see Figure 5). Each area of the plot corresponds to different types 
of X-ray sources. Figure 5 further demonstrates the strange behavior of this source, as it 
has detections in very different locations on the graph. 



4. Optical Observations and Analysis 



The location of rl-25 has been observed with the ACS camera onboard HST (see Table 
10). In both cases, images are only available in the F435W filter (approximately equal to 
the B filter in ground based systems). We have 4 observations aquired in the period 2004- 
2007, so the data are useful for looking for transient or variable sources. Only one Chandra 
observation (ObsID 8183) was taken within a week of an ACS observation (j9ju06010) of the 
source. However, the source was not detected in ObsID 8183. 

All data were obtained using th e Wide Field Channel (WFC), which has a 202" x 202" 
field of view jMavbhate et al.ll201oh. _ Each observation was carried out in the standard four 
pointing dither pattern (IMaybhate et al.ll2010[ ). These individual images were then combined 
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using the PyRAF task MultiDrizzle (IFruchter et al.ll2009[ ). which also removes cosmic rays 
and corrects the geometrical distortion which results from the orientation of the ACS with 
respect to the HST focal plane. We chose not to apply an automatic background subtraction 
in order to perform photometry, since there is a steep gradient in the diffuse light this close to 
the center of M31 making background estimation unreliable. Finally, we utilized the ability 
of MultiDrizzle to resample the spatial scale of the image, resulting in a final pixel scale of 
0.025'7pix. 

The World Coordinate System (WC S) in HST images can be offset from standard refer- 
ence frames by as much as one arcsecond ( Maybhate et al.lboiol ). To improve the astrometry, 
we regis tered the final drizzl ed images to the WCS of the Local Group Survey (LGS) images 
of M31 (IMassey et al.l 120061 ). Stars common to both images were identified, and their cen- 
troid positions calculated using the IRAF task imcentroid. We then used the task ccmap in 
IRAF to update the WCS of the HST images. The final rms (Icr) errors in the alignment 
were of order 0.006" in RA, and 0.002" in declination. We note that the RMS errors on the 
alignment of the HST images to the WCS of the LGS survey were always smaller than 0.01" 
(which is smaller than one pixel in the rescaled images). 

We also aligned the deepest Chandra observation (ObsID 1575) with the WCS of the 
LGS images using the same procedure applied to the HST images. We found an alignment 
error of 0.109" in RA and 0.149" in DEC. The centroid position in the corrected WCS for the 
source is RA:00:42:47.90, DEC:+41: 15:49.99, with errors of 0.08" (RA) and 0.07" (DEC). 
To get the final positional error, we added the alignment and centroid errors in quadrature. 
The final RA error is 0.13" and final DEC error is 0.17" (these are la errors). 

Since the WFC images cover such a large area and contain so many stars, we extracted 
subimages 1000 pixels on a side centered on our source of interest and performed pho tometry 
on those. We used the DAOPHOT II and ALLSTAR routines Jstetson et al.lll99oh to find 
and photometer stars in the images. Stars suitable for calculating a point spread functions 
(PSFs) were identified by hand to avoid problems due to crowding. The final measured 
counts for each star were converted first to count rate by dividing by the exposure time of 
each observation, and then to AB system m agnitudes using the conversion factors in the 
ACS users handbook jMavbhate et al.lboioh . 



We show the final reduced image for the source in Figure 6, superimposed with the 
X-ray 3cr positional error ellipse. The ellipse is drawn using the 3(T RA error (0.39") as the 
semi-major axis and the 3(J DEC error (0.51") as the semi-minor axis. The location of rl-25 
is extremely crowded, and being so close to the core is subject to a high background of diffuse 
light. In our series of four images, two (j9ju01010 and j9ju06010) are also significantly deeper, 
which affects the completeness of the stars we can detect. The small size of the Chandra 
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enor circle does however simplify the analysis of the photometry, since very few sources are 
inside the X-ray 3cr error ellipse. Examining the images, a number of sources are detected 
inside the error circle, although most are unresolved. We note that there are no catalog 
stars within the Chandra error ellipse. In fact, only one source is resolved inside the error 
ellipse in all four images by the photometry source detection algorithm. We have marked 
this source with a red circle in Figure 6. 

The single resolved source in Figure 6 has observed magnitudes of 24.45 ± 0.06, 24.27 ± 
0.05, 24.31 ± 0.04 and 24.41 ± 0.04 in each of the four images, whe re the quoted error i s the 



la magnitude error. Using the standard E{B — V) = 0.062 for M31 (jSchlegel et al.lll998l ). the 



extinction in the F435W filter is 0.785 magnitudes. Thus, the star's Mr is between -0.80 and 



-0.98 in the four observations, assuming (m— M)o = 24.47 (lHollandlll998l ). This demonstrates 
that, within the uncertainties, there is no evidence of variability in this object. Although 
other sources are picked up by the DAOphot detection algorithm inside the error ellipse in 
some images, these additional detections can be accounted for by the longer exposure time, 
or are unreliable due to crowding. 



(IGrupe et al.ll2010l ) looked at the spectral energy distribution of 92 active galactic nuclei 
that had soft X-ray spectra. The AGN they studied had comparable X-ray count rates to 
rl-25. However, the AGN were much brighter in the B band (14 < rriB < 18 ) than any 
source in Figure 6. For this reason, we are confident that rl-25 is not an AGN. 

With no color information, we cannot determine what the object marked in Figure 6 is 
for certain. If it is a star, it would correspond to a late B type with bolometric luminosity 
of 10^^ ergs s~^. We note that it is too luminous in the B band to be a red giant, and too 
dim to be a red supergiant. 



5. Models 

The state changes observed in rl-25 are extremely unusual for a SSS, and quite different 
than state changes observed in other classes of X-ray source. In this section, we consider 
a number of physical models with the goal of uncovering the nature of the X-ray emitting 
source. In order for any model of rl-25 to be successful, it must be able to explain all of the 
observed features of the source. These features include the source's appearance as an SSS- 
HR source in observation 1575 with ACIS-S, with kT ~ 130 eV and a 0.3-8 keV luminosity 
of 4 X 10^^ erg s~^. We wish to emphasize once again that even though rl-25 had an effective 
temperature in excess of 100 eV when detected as an SSS (much higher th an most SSSs), 



it satisfied the strictest SSS criterion as defined in iDi Stefano et al.l (120041 ). That is, the 
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detection in observation 1575 (when the source was an SSS) had no hard counts, medium 
counts consistent with zero, and a 3cr detection in the soft band. 



Furthermore, the model must explain the subsequent detections of the source as a ~250 
eV source, with higher luminosity than in the soft state (~10^'' erg s^"*^). Lastly, the model 
must be consistent with an optical counterpart with F435W magnitude fainter than ~24.3 
{Mb fainter than ~ -1). We note that white dwarf, neutron star, and black hole SSSs 



are consistent with this optical restraint. CAL 87, with Mb ~ 0.5 (e.g., iGreinei 



Deutsch et al.lll996l ). is thought to be a stellar- mass black hole system fe.g... ICowlev et al. 



2000 



2007 



1990l ). RX J0720.4-3125, an isolated neutron star SSS, has a Ms ~ 18 JXaplan et al. 
references therein). RX J0439. 8-6809, a white dwarf SSS, has a Mb ~ 2.7, which is also 
consistent with the optical constraint on rl-25. 

The following subsections outline white dwarf, neutron star and black hole models, and 
compare their features to the observed properties of rl-25. 



5.1. White Dwarfs 



White dwarfs that have recently experienced novae have temperatures and luminosities 
that can make them detectable as SSSsH Many SSSs detected in M31 were recent novae 
(e.g., Pietsch et al. 2005 , 2007 : Stiele et al.lbolO , references therein). When the white dwarf 
cools, some novae ar e detected as har der X-ray sources, but the X-ray luminosity is around 
Iq32 _ j^gss gj^gg g-i ( jSala et al.l 120101 ). In contrast, n ovae in a supersof t state are detected 
with A;T ~ 50 eV and Ly of 10^^ - lO^^ ergs s'^ fe.g.. Istiele et al.l b 



Stiele et al.ll2010l . references therein). 



Thus rl-25, even in its softest state (with kT ~130 eV and Lx of 4 x 10^^ erg s~^), is too 
hard to be consistent with the very soft emission detected in typical supersoft novae. Also, 
the source is not consistent with the harder states of novae. That is rl-25, even in its hardest 
state (~250 eV and Lx of 1.1 x 10^^ erg s~^) is softer and more luminous than novae in 
their hardest states. Moreover, the harder states of rl-25 are more luminous than the softer 
ones, in contradiction to what happens in novae. We therefore turn to models in which a 
white dwarf accretes matter at high rates. 

When a white dwarf accretes mass at a high enough rate that the incoming matter can 
experience nuclear burning, the white dwarf can appear as an SSS. In this case, the source 



specific post-nova system will be detectable as an SSS only if the white dwarf stays hot enough 
to be emitting as an SSS until after the optical depth has decreased enough to let radiation escape (e.g., 
Sala fc Hernaii3l2005h . 
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will not be in a hard state. The copious energy we receive from such sources is provided by 
nuclear burning, rather than accretion. The effective radii are comparable to the white dwarf 
radii, so the emission can be characterized by values of fcT in the range of tens of eV for 
low-mass white dwarfs, and ~ 100 eV for white dwarfs approaching the Chandrasekhar mass 
(Mc). For each w hite dwarf mass, nuclear burning can occur only within a narrow range 



of accretion rates ( iDi Stefandl2010l ). These rates are very high: ~ 10 ''Mcd for a solar- mass 



white dwarf, ten times higher for a white dwarf with mass near Mc ( iDi Stefandl2010l ). At 
such rates of infall, accretion alone produces luminosities in the range of 10^^ — 10^^ ergs s~^, 
typically a few percent of the total energy of the system. 



Consider a case, when its accretion rate places a white dwarf near the lower end of 
the steady-burning region, or just below it. In this case, nuclear burning may be episodic. 
During and just after nuclear-burning episodes, the emission is dominated by soft emission. 
As the white dwarf cools, however, it becomes less luminous and the emission is dominated 
by accretion. Although at high rates of accretion, the emission is expected to be softer than 
typical for low-accretion-rate white dwarfs, such as cataclysmic variables, it can nevertheles s 



be harder than typical of SSSs (e.g. JPopham fc Narayanlll995l : [Patterson fc Raymondlll985l ). 



Thus, the source could appear to be quasisoft or hard. If the donor is a giant or a Roche- 
lobe filling star in a circular orbit, the accretion rate should could continue to be high. The 
source will continue to be detected as a harder source. Nuclear burning episodes occurring 
at intervals ranging from months to decades would make the source more luminous and 
detectable as an SSS. 

The points made above about the quasi-steady nuclear burning white dwarf model are 
illustrated in Figure 7. The figure is a plot of kT vs. the logarithm of bolometric luminosity, 
LOG[L] for various quasi-steady nuclear burning white dwarfs (green points) along with the 
rl-25 spectra (black, blue, and red points). The solid and open black points represent the 
XSPEC spectral fits of rl-25. The solid black points represent the fits where we assume 
Lx = L. The X-ray luminosity, however, is not equivalent to the bolometric luminosity. We 
assume that the X-ray luminosity represents, at a minimum, a quarter of the bolometric 
luminosity. Thus, we plot the open black circles which represent XSPEC fits assuming 
Lx = 0.25L. The solid blue points, which represent the spectral fits for observations with 
fewer counts (spectra estimated using the pimms method described in §3.2 ), are included 
for completeness. They represent the pimms spectral estimates assuming Lx = L. The red 
points represent the pimms results assuming Lx = 0.25L. However, as noted in §3.2, these 
spectral fits are estimates. For this reason, we point out that the XSPEC results are more 
reliable. In fact, the error bars are so large for the pimms results that they cover most of 
the plot. The pimms uncertainties are excluded for this reason. The plot clearly indicates 
that the quasi-steady nuclear burning white dwarf model does not fit rl-25, as the data are 
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to hard and/or dim to fit the model. That is, no XSPEC points and only one of the pimms 
points fall in the range of quasi-steady nuclear burning. 



5.2. Neutron Stars 

Isolated neutron stars and neutron stars in quiescent low-mass X-ray binaries (qLMXBs) 
have been observed with spectra in the SSS or QSS range, but they are typically 3-5 orders of 



magn it ude less luminous t h an th e "classical" SSSs and QSSs (e.g.. iHaberll 120071 : iPires et al. 
20091 ). iKylafis &: Xilourid ( 119931 ) showed, however, that the combination of gravitational 
and magnetic effects can cause neutron stars accreting near the Eddington limit to produce 
SSS-hke spectra. In such cases, the luminosities and temperatures would give the accreting 
neutron stars the appearance of "classical" SSSs. Nevertheless, neutron-star models for SSSs 
have received relatively little attention, perhaps because of the success of white dwarf models, 
coupled with what appeared to be the fine tuning required for a neutron star to achieve the 
appropriate effective radius. 

It is therefore worth noting tha t there i s obser vational support for the idea that neutron 
stars can produce very soft spectra. iHughed (Il994l ) discovered a transient pulsar in the Small 
Magellanic Cloud that has an unpulsed, highly luminous (near Eddington) soft ~ 35 eV) 
component. If such a system is viewed from an angle at which the hard radiation is not 
detected, it would have t he properties asso ciated with the "classical" SSSs first discovered 
in the Magellanic Clouds ( jLong et al.lll98ll ). 



Here we point out that the discovery of QSSs provides reasons to revisit neutron star 
models. Consider a neutron star producing near-Eddington soft emission. If this emission 
emanates from a nearly spherical photosphere with radius equal to the Alfven radius, then 



A; T = 104 eV 



10"G 

5., 



10 km 



R 



NS 



M 

Mo 



L 



10^''ergss 



-1 



(4) 



In this expression, R^s is the radius of the neutron star, Bg is the value of magnetic field 
on the surface, and M and L represent the neutron star's mass and luminosity, respectively. 

An interesting feature of this expression is that there are ranges of reasonable values 
of the physical parameters in which the /cT is in the range expected for SSSs or QSSs. 
Furthermore, for a specific neutron star, the effective temperature depends on L which can 
change as the accretion rate changes. Since changes in accretion rate are common, we may 
therefore expect the effective temperatures of some neutron stars to change. Depending on 
the physical parameters, these changes could produce transitions from SSS to QSS states. 
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Another set of factors support neutron star models for some QSSs and state-changing 
QSSs. This is that, if there is luminous soft emission, there must also be a hard component 
generated by matter accreting, which is channeled onto the poles by the neutron star's strong 
magnetic field. The hard emission is likely to be beamed. This model therefore predicts that 
a large population of accreting neutron stars will contain some with a pulsed hard component. 
If this component is not the dominant one, then the neutron star may appear as a QSS, even 
if emission from the Alfven sphere is SSS. The diagnostic for the neutron star model would 
be periodicity in the arrival times of the harder photons. 

Figure 8 shows two plots of kT vs. LOG[L]. The top pannel shows several curves for 
various neutron star models. Each curve differs by a factor of ten in Bs- The bottom plot is 
of the rl-25 spectra along with the Bg = 10^" G curve. The black (XSPEC) and blue (pimms) 
points were plotted in the same manner as in Figure 7 (see §5.1), except that we assumed 
for neutron star models Lx ~ L only. The distribution of points is roughly consistent with 
what is expected for this neutron star model. The XSPEC points seem to follow the same 
trend as the curve, with some variation. There is more variability in the pimms data, but 
the uncertainties are large for these spectral estimates. 

It is important to note that if we can identify QSSs and state changers in the Milky 
Way or in the Magellanic Clouds, then neutron star models are testable. The reason is that 
neutron stars produce pulsed radiation and/or hard bursts. This is how we have verified 
the neutron-star n atures of LMXBs in the Galaxy's gl obular clusters. Both bursts (e.g.. 



Lewin et al.l Il993l ) and pulsed (e.g., IZhelezniakovl Il98l[ ) radiation would be detectable in 



nearby neutron stars, even if the radiation is dominated by an unpulsed soft component. 



5.3. Black Holes 



Accreting black holes can exhibit thermal-dominant states in which the emission is 
dominated by a component em anating from the inner portion of the accretion disk (e.g., 
Remillard fc McClintockll2006[ ). SSSs and QSSs have therefore both been suggested as pos- 
sible black holes. In fact, the most well-know state changer is MlOl-ULX-1, an ultralurninous 
SSS that has been detec t ed also in high Q SS and low-hard states (e.g., iKong et al.l 12004 : 
Kong fc Pi Stefanoll2005l : iMukai et al]l2005h . MlOl-ULX-1 is almost certainly a black hole. 
Its mass could be either in the range typical of Galactic stellar-mass black holes or else in 
the higher range (50 — 10^ Mq) suggested for intermediate-mass black holes. 

The luminosity of rl-25 is 1-3 orders of magnitude smaller than the luminosities mea- 
sured for MlOl-ULS-1 when it is in a soft state. It is therefore highly unlikely to be an 
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intermediate-mass black hole. In fact, if the luminosity is less than roughly a percent of 
the Eddington luminosity, then the inner disk will not be optically thick and the emission 
will not be thermal. This suggests that, if this source is a black hole, it is more h kely to 



be of stellar mass. The top panel of Figure 9, first shown in lDi Stefano et al.l ( 120101 ). shows 
that QSS emission is expected in the thermal-dominant state of black holes with mass below 
~ 100 Mq. The radius of the inner disk would determine the value of the effective temper- 
atures; the spectrum could be either QSS or SSS. At lower rates of accretion, the emission 
would be hard. 

The bottom panel of Figure 9 shows the rl-25 spectra along with the lOM© and lOOM© 
curves. The black (XSPEC) and blue (pimms) points were plotted in the same manner as 
in Figure 8 (see §5.2). The distribution of the points with XSPEC fits is roughly consistent 
with what is expected for a black hole of roughly lOM©. The distribution of points with 
fewer counts (pimms fits) is broader, but the uncertainties are larger. 



Conclusion 



We have tracked the long-term behavior of the M31 X-ray source rl-25. First observed 
by ROSAT on 1990-07-24, then by both Chandra and XMM-Newtona.nd most recently by 
Swift, rl-25 is one of the best-studies soft X-ray sources. There are 86 public Chandra 
observations of the source through June 2, 2009, with 45 detections . For XMM- Newton, 
there are 26 public observations of the source, with IStiele et al.l (120101 ) reporting detections 
in only the 2004 data. The detections of the source start in 1999 and continue through 2009. 

By doing this we have documented the fact that rl-25 has transitioned from an SSS to 
a harder, QSS state. In the SSS state its estimated X-ray luminosity is a few times 10^^ 
ergs s~^, and the luminosity appears to be higher, but not much over 10^^ ergs in the 
harde r state. Only oiie other X-ray source, MlOl-ULX-1, has wel l-studied state changes 
(e.g., iKong et al.l 120041 : iKong fc Pi Stefanol I2OO5I iMukai et al.ll2005h . While MlOl-ULX-1, 



which has been observed with X-ray luminosity as high as lO^'' — 10^^ ergs s~^, is almost 
certainly a black hole, the nature of rl-25 is more difficult to establish, because its luminosity 
range is consistent with white dwarf, neutron star, or black hole accretors. 

Whatever its nature it is different from anything we have observed. We have shown 
that the observed behavior is consistent with a black hole accretor with a mass in the IOMq 
range. In this case, our observations of rl-25 have all found it to be in a thermal- dominant 
state. The inner disk radius would have been larger in the SSS state. If rl-25 is a black hole 
with a mass of approximately IOMq, it could be more luminous in future observations, if 
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the donor star is able to contribute mass at a higher rate. Should the luminosity approach 
the Eddington luminosity, the system would be unlikely to remain in the thermal dominant 
state, and hard emission could be detected. Similarly, if the luminosity falls below ~ 1% of 
the Eddington value, the spectrum would likely be hard. 

We have also shown that the observed behavior of rl-25 is consistent with a neutron star 
accretor with Bg = 10^^ G. We assume that the magnetic field should be constant over the 
short interval of observations of the source. The model suggests that the harder states are 
more luminous than the softer ones, which is consistent with the rl-25 spectra. We note that 
neutron star models are testable if we can find state changers and QSSs in the Milky Way 
or Magellanic Clouds, as both bursts and pulsed radiation would be detectable in nearby 
neutron stars.. 

White dwarf models are the least likely fit for rl-25. The source does not seem to exhibit 
behavior of a post-nova system, as its spectrum is harder than novae that are SSSs (even 
when rl-25 is in its soft state). The source is also too luminous and soft to be consistent 
with novae in there harder states. Furthermore, quasi-steady nuclear burning white dwarf 
models do not fit the data. We have shown that quasi-steady nuclear burning models can 
produce both supersoft and quasisoft radiation, but the model is too soft and/or luminous 
to fit the rl-25 spectra. 

We note that there are other state-changing sources in external galaxiesH For example, 
there are nine state changing sources in nearby spiral galaxy M33. If we study a large enough 
sample of state-changers, we are likely to find examples of all three (white dwarf, neutron 
star, and black hole) models. Continued monitoring of these sources will play an important 
role in testing these models. It is also important to identify QSSs and state-changers in 
the Magellanic Clouds and Galaxy, where many test of the nature of the accretors can be 
conducted. 

We would like to thank the Swift team for approving our ToO request (Target ID: 
35336). This research has made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, Caltech, under contract with the National 
Aeronautics and Space Administration. This research was supported by HST Grant AR- 
10948. 01-A-O and the Smithsonian Institution IR & D award. 
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Fig. 1. — : The figure shows the location of rl-25 (the green circle) in M31 in different 
observations of rl-25. Top Left: Merged HRC data of the source. Top Right: ACIS-S ObsID 
1575 of the source. Bottom: ACIS-I ObsID 4721 of the source. The location of the source is 
the same in the two ACIS observations. Although the HRC data do show an offset of ~0.5" 
from the rl-25 ACIS position, other sources in the field do as well. Thus, the offset comes 
from a shift in the coordinate frame. 
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Fig. 2. — : Source (solid) and background (dashed) apertures for rl-25, displayed on a ~432 
kilosecond merged Chandra HRC-I observation of M31. 
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Fig. 3. — : Light curve for rl-25. Bottom Panel: Count Rate in the Soft (0.1-1.1 keV), 
Medium (1.1-2 keV), and Hard (2-7 keV) bands vs. Obsdate. Middle Panel: Total counts 
for ACIS observations of rl-25 vs. Obsdate. Bottom Panel: TCR in ACIS-I units (in the 
energy range of 0.1-7 keV) vs. Obsdate. See text for how we converted ACIS-S and HRC-I 
observations to ACIS-I units. 
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Fig. 4. — : rl-25 spectral analysis. The blue and black symbols represent spectral results 
from pimms and XSPEC, respectively. Plots assume Nh = 1 x 10^^ cm~^. Top Panel: 
Plot of logarithm of X-ray luminosity, LOG[Lx], vs. Obsdate for ACIS observations that 
detected rl-25 (Lx in the 0.1-7 keV range). Middle Panel: Plot of kT vs. Obsdate for ACIS 
observations that detected rl-25. Bottom Panel: Plot of LOG[Lx] vs. kT. 
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Fig. 5. — : Only detections of rl-25 with 25 or more counts were included on this color-color 
graph. The numbers refer to the observations in Table 9. The ellipses were taken from 



Prestwich et al.l (120031 ). Ellipse A is generally associated with SNRs, Ellipse B with LMXBs, 
Ellipse C with absorbed sources, and Ellipse D with HMXBs. The bands used for this plot 
were: S,0.3-l keV; M,l-2 keV; H,2-8 keV; and T,0.3-8 keV. 
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Fig. 6. — : HST ACS images of the location of rl-25, obtained in the F435W filter on four 
separate dates (see Table 13). Top Left: Observation j8vp03010, Top Right: observation 
jSvpOSOlO, Bottom Left: Observation jQjuOlOlO, Bottom Right: Observation j9ju06010. Each 
image is 3" on a side. The black ellipse denotes the 3a X-ray positional error ellipse of rl-25. 
The ellipse has a semi-major axis {3a RA error) = 0.39" and a semi-minor axis {3a DEC 
error) = 0.51 ". The red circle shows the only consistently detected, resolved source found 
by DAOphot. 
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Fig. 7. — : Plot of kT vs. the logarithm of the bolometric luminosity, LOG[L]. The black 
and blue points are the average kT and LOG[L] for XSPEC and pimms spectral fits of rl-25, 
respectively (assuming = 1 x 10^^ cm~^). The solid black points represent XSPEC fits 
assuming Lx = L, and the open black points assume Lx = 0.25L for XSPEC fits. The blue 
points represent the pimms results assuming Lx = L. The red points correspond to pimms 
results assuming Lx = Lx = 0.25L (see text). The error bars represent the range of kT 
and Lx that fit each observation with XSPEC results (see Table 8). The error bars are not 
plotted for the pimms points as the uncertainties are too large (see text). The green points 
represent the kT and LOG[L] for various quasi-steady nuclear burning white dwarfs. 
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Fig. 8. — : Plot of kT versus L for the neutron star model described by Equation 4. Top 
Panel: Each curve corresponds to a fixed value of the magnetic field on the surface of the 
neutron star; Bg changes by a factor of ten between curves. Since the field is not expected 
to vary this much over short times, the system should evolve along the curves. Thus, if the 
luminosity increases or decreases, so does the temperature. Bottom Panel: The Bg = 10^° G 
curve is plotted with the rl-25 spectra. The black and blue points are the average k T and 
LOG[L] for XSPEC and pimms spectral fits of rl-25, respectively (assuming Nh = 1 x 10^^ 
cm~^). The error bars represent the range of kT and Lx that fit each observation (see Table 
8). The error bars are not plotted for the pimms points as the uncertainties are too large 
(see text). 
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Fig. 9. — : Figure taken from iDi Stefano et al.l (120101 ) . It is a plot of kT versus LOG[L] 
for the inner portion of the accretion disk around black holes. Top Panel: Each pair of 
two curves of a single color corresponds to a fixed black hole mass which labels the regions 
between the curves. The upper curve of each color corresponds to a disk with inner radius 
GMbhG/c^, while the lower curve corresponds to an inner disk with 3 times the radius. The 
point at the bottom (top) of each curve corresponds to the luminosity of the inner disk being 
l%LEddington {^0%LEddington) ■ BottoTTi Panel: The lOM© and IOOM0 curves are plotted with 
the rl-25 spectra. The black and blue points are the average kT and LOG[L] for XSPEC 
and pimms spectral fits of rl-25, respectively (assuming Nh = 1 x 10^^ cm^^). The error 
bars represent the range oi kT and Lx that fit each observation (see Table 8). The error 
bars are not plotted for the pimms points as the uncertainties are too large (see text). 



-30- 



Table 1. Chandra Observations of r 1-25 







Exposure Time 


Observation 




Detected 


ObsID 


Instrument 


seconds 


Date 


OBSMJD 


rl-25? 


303 


ACIS-I 


11862 


1999-10-13 


51464.2 


No 


267 


HRC-I 


1270 


1999-11-30 


51512.8 


No 


305 


ACIS-I 


4129 


1999-12-11 


51523.2 


No 


268 


HRC-I 


5178 


1999-12-23 


51535.1 


No 


306 


ACIS-I 


4132 


1999-12-27 


51539.7 


No 


269 


HRC-I 


1206 


2000-01-19 


51562.9 


No 


307 


ACIS-I 


4113 


2000-01-29 


51572.6 


No 


270 


HRC-I 


1470 


2000-02-13 


51587.2 


No 


308 


ACIS-I 


4012 


2000-02-16 


51590.1 


No 


271 


HRC-I 


2464 


2000-03-08 


51611.8 


No 


272 


HRC-I 


1203 


2000-05-26 


51690.7 


No 


309 


ACIS-S 


5089 


2000-06-01 


51696.1 


No 


273 


HRC-I 


1191 


2000-06-21 


51716.2 


No 


310 


ACIS-S 


5074 


2000-07-02 


51727.9 


No 


311 


ACIS-I 


4894 


2000-07-29 


51754.0 


No 


275 


HRC-I 


1188 


2000-08-18 


51774.4 


No 


312 


ACIS-I 


4666 


2000-08-27 


51783.8 


Yes 


276 


HRC-I 


1183 


2000-09-11 


51798.9 


No 


277 


HRC-I 


1187 


2000-10-12 


51829.2 


No 


278 


HRC-I 


1182 


2000-11-17 


51865.2 


No 


1581 


ACIS-I 


4404 


2000-12-13 


51891.1 


Yes 


1854 


ACIS-S 


4694 


2001-01-13 


51922.4 


Yes 


1569 


HRC-I 


1169 


2001-02-01 


51941.0 


No 


1582 


ACIS-I 


4303 


2001-02-18 


51958.7 


No 


1583 


ACIS-I 


4903 


2001-06-10 


52070.8 


Yes 


1570 


HRC-I 


1184 


2001-06-10 


52070.9 


No 


1577 


ACIS-I 


4916 


2001-08-31 


52152.5 


No 


1575 


ACIS-S 


37664 


2001-10-05 


52187.0 


Yes 


1912 


HRC-I 


46732 


2001-10-31 


52214.0 


Yes 


1585 


ACIS-I 


4884 


2001-11-19 


52232.8 


No 


2904 


HRC-I 


1182 


2001-11-19 


52232.9 


No 


2895 


ACIS-I 


4881 


2001-12-07 


52250.8 


No 


2897 


ACIS-I 


4900 


2002-01-08 


52282.1 


No 


2905 


HRC-I 


1094 


2002-01-16 


52290.2 


No 


2896 


ACIS-I 


4903 


2002-02-06 


52311.9 


No 


2898 


ACIS-I 


4900 


2002-06-02 


52427.9 


No 


2906 


HRC-I 


1188 


2002-06-02 


52427.9 


No 


4360 


ACIS-I 


4906 


2002-08-11 


52497.7 


No 


4678 


ACIS-I 


3894 


2003-11-09 


52952.3 


Yes 


4679 


ACIS-I 


3820 


2003-11-26 


52969.9 


Yes 


4680 


ACIS-I 


4198 


2003-12-27 


53000.3 


Yes 


4681 


ACIS-I 


4094 


2004-01-31 


53035.8 


No 


4682 


ACIS-I 


3945 


2004-05-23 


53148.7 


Yes 


4719 


ACIS-I 


4123 


2004-07-17 


53203.9 


Yes 


4720 


ACIS-I 


4108 


2004-09-02 


53250.6 


Yes 
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Table 1 — Continued 







Exposure Time 


Observation 




Detected 


ObsID 


Instrument 


seconds 


Date 


OBSMJD 


rl-25? 


4721 


ACIS-I 


4132 


900A 1 Ozl 


KQOQO Q 
OOZoZ.y 


leo 


4722 


ACIS-I 


3894 


900zL 1 '^l 


'^'^'^OQ 1 
uoouy . -L 


Yes 


4723 


ACIS-I 


4035 


900zL 1 9 0'^ 


53344.4 


Yes 


5925 


HRC-I 


46311 


9004-1 9-.0R 


uootu .o 


Yes 




HRC-I 


zuuoo 


9004 1 9 97 


DOOUIJ.O 


leo 


5926 


HRC-I 




9004-1 9-97 


uoouu.o 


Yes 


6202 


HRC-I 




900"=* 01 9R 




leb 


5927 


HRC-I 


27000 


900^^ 01 9^^ 


53398.8 


Yes 


5928 


HRC-I 


44856 


900^^ 09 91 




Yes 


7136 


ACIS-I 


3991 


900fi 01 OR 


53741.8 


Yes 


7137 


ACIS-I 


3952 


900R-0'^-9R 


'^'^KSI 9 

UOOO J- 


Yes 


7283 


HRC-I 


19942 


900R-0R-0'=1 


tjooy J. .o 


Yes 


71 "^Si 

/ -LOO 


A cm T 


A^ 07 


900R OR OQ 




leb 


7139 


ACIS-I 


3987 


900R 07 "^l 




Yes 


7140 


ACIS-I 


4117 


900R OQ 94 


•^4009 ^ 


Yes 


7284 


HRC-I 




900R-0Q-'^0 
^uuu uy ou 


•^400^ Q 


Yes 


7285 


HRC-I 


18517 


900R- 11 - m 


^40^9 ^ 


Yes 






OQOQQ 


900R 1 9 04 


■^407*^ Q 
u^u 1 o.y 


leo 


8183 


ACIS-I 


3965 


9007-01-14 


54114.9 


No 


8184 


ACIS-I 


4140 


9007 09 1 4 

ZiVjVJ I UZ -L'± 


54145.1 


No 


8185 


ACIS-I 


3965 


9007 0*^ 1 

ZrUU I -UO -LU 


•^41 RQ 9 


No 


7286 


HRC-I 


18708 


9007-0'^- 1 1 


54170.6 


No 


7068 


ACIS-I 


7693 


9007-0R-09 




Yes 


8191 


ACIS-I 


3959 


2007-06-18 


54269.3 


No 


8192 


ACIS-I 


4073 


2007-07-05 


54286.5 


Yes 


8193 


ACIS-I 


4129 


2007-07-31 


54312.1 


Yes 


8194 


ACIS-I 


4033 


2007-08-28 


54340.5 


Yes 


8195 


ACIS-I 


3965 


2007-09-26 


54369.6 


Yes 


8186 


ACIS-I 


4134 


2007-11-03 


54407.2 


Yes 


8526 


HRC-I 


19944 


2007-11-07 


54411.6 


Yes 


8527 


HRC-I 


19981 


2007-11-17 


54421.8 


Yes 


8187 


ACIS-I 


3837 


2007-11-27 


54431.2 


Yes 


8528 


HRC-I 


19975 


2007-11-28 


54432.8 


Yes 


8529 


HRC-I 


18923 


2007-12-07 


54441.6 


Yes 


8530 


HRC-I 


19915 


2007-12-17 


54451.5 


Yes 


9529 


ACIS-I 


4113 


2008-05-31 


54617.5 


Yes 


9520 


ACIS-I 


3962 


2008-12-31 


54463.7 


Yes 


10551 


ACIS-I 


3965 


2009-01-09 


54841.0 


No 


10552 


ACIS-I 


3962 


2009-02-07 


54869.4 


No 


10553 


ACIS-I 


4104 


2009-03-11 


54901.6 


Yes 


10554 


ACIS-I 


4036 


2009-05-29 


54980.1 


No 
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Table 2. Photometry for rl-25 



Soft Medium Hard Total SCR MCR HCR TCR ACIS-I Rate 

OBSM.TD Counts Counts Counts Counts 10~3s-l IQ-^s-^ 10~3s-l IQ-^s-l lO-^s-^ 



51464 


2 




















< 1.63 


< 1.63 


51512 


8 




















< 4.66 


< 5.87 


51523 


2 




















< 2.70 


< 2.70 


51535 


1 




















< 1.85 


< 2.33 


51539 


7 




















< 2.70 


< 2.70 


51562 


9 




















< 4.90 


< 6.17 


51572 


6 




















< 1.44 


< 1.44 


51587 


2 




















< 5.29 


< 6.66 


51590 


1 




















< 1.47 


< 1.47 


51611 


8 




















< 3.73 


< 4.70 


51690 


7 




















< 6.46 


< 8.13 


51696 


1 




















< 4.20 


< 2.17 


51716 


2 




















< 4.97 


< 6.26 


51727 


9 




















< 2.71 


< 1.40 


51754 























< 2.16 


< 2.16 


51774 


4 




















< 6.53 


< 8.22 


51783 


8 


6.92 ± 3.80 


0.00 ± 1.87 


0.00 ± 


1.87 


6.92 ± 3.79 


1.50 ± 


0.82 


0.00 ± 0.41 


0.00 ± 0.41 


1.60 ± 0.82 


1.50 ± 0.82 


51798 


9 




















< 6.76 


< 8.51 


51829 


2 




















< 4.98 


< 6.27 


51865 


2 




















< 9.57 


< 12.05 


51891 


1 


5.83 ± 3.60 


1.00 ± 2.32 


1.00 ± 


2.32 


7.83 ± 3.96 


1.34 ± 


0.83 


0.23 ± 0.53 


0.23 ± 0.53 


1.80 ± 0.91 


1.80 ± 0.91 


51922 


4 


19.75 ± 5.56 


3.00 ± 2.94 


0.00 ± 


1.91 


22.66 ± 5.88 


4.26 ± 


1.20 


0.65 ± 0.63 


0.00 ± 0.41 


4.89 ± 1.27 


0.79-1.48 


51941 























< 7.97 


< 9.36 


51958 


7 




















< 1.79 


< 1.79 


52070 


8 


6.75 ± 3.79 


0.00 ± 1.87 


0.00 ± 


1.87 


6.75 ± 3.79 


1.39 ± 


0.78 


0.00 ± 0.38 


0.00 ± 0.38 


1.39 ± 0.78 


1.39 ± 0.78 


52070 


9 




















< 9.56 


< 11.22 


52152 


5 




















< 2.45 


< 2.45 


52187 





169.67 ± 14.16 


14.33 ± 4.98 


0.00 ± 


1.88 


183.76 ± 14.71 


4.55 ± 


0.38 


0.38 ± 0.13 


0.00 ± 0.05 


4.92 ± 0.39 


1.71-2.14 


52214 























0.46 ± 0.13 


0.147-0.686 


52232 


8 




















< 1.86 


< 1.86 


52232 


9 




















< 6.50 


< 11.76 


52250 


8 




















< 1.61 


< 1.61 


52282 


1 




















< 1.59 


< 1.59 


52290 


2 




















< 7.30 


< 13.21 


52311 


9 




















<1.59 


<1.59 


52427 


9 




















< 1.21 


< 1.21 


52427 


9 




















< 8.19 


< 14.82 


52497 


7 




















< 1.58 


< 1.68 


52952 


3 


7.92 ± 3.96 


2.00 ± 2.66 


0.00 ± 


1.87 


9.92 ± 4.28 


2.06 ± 


1.03 


0.52 ± 0.69 


0.00 ± 0.49 


2.58 ± 1.11 


2.58 ± 1.11 


52969 


9 


6.67 ± 3.79 


().()() ± 1.87 


2.00 ± 


2.66 


8.67 ± 4.13 


1.77 ± 


1.01 


0.00 ± 0.50 


0.53 ± 0.71 


2.30 ± 1.10 


2.30 ± 1.10 


53000 


3 


8.75 ± 4.13 


3.00 ± 2.94 


1.00 ± 


2.32 


12.75 ± 4.71 


2.11 ± 


0.99 


0.72 ± 0.71 


0.24 ± 0.56 


3.07 ± 1.13 


3.07 ± 1.13 


53035 


8 




















< 3.22 


< 3.22 


53148 


7 


12.92 ± 4.71 


10.00 ± 4.28 


0.00 ± 


1.87 


22.92 ± 5.87 


13.31 ± 


4.85 


10.30 ± 4.41 


0.00 ± 1.93 


23.62 ± 6.05 


23.62 ± 6.06 


53203 


9 


18.67 ± 6.46 


15.92 ± 6.09 


0.00 ± 


1.87 


34.68 ± 6.98 


4.63 ± 


1.36 


3.95 ± 1.26 


0.00 ± 0.46 


8.68 ± 1.73 


8.58 ± 1.73 


53250 


6 


24.92 ± 6.08 


34.00 ± 6.89 


2.92 ± 


2.94 


61.83 ± 8.92 


6.10 ± 


1.49 


8.33 ± 1.69 


0.72 ± 0.72 


15.14 ± 2.18 


16.14 ± 2.18 


53282 


9 


29.00 ± 6.46 


33.83 ± 6.90 


2.92 ± 


2.94 


66.75 ± 9.17 


7.11 ± 


1.68 


8.30 ± 1.69 


0.72 ± 0.72 


16.12 ± 2.25 


16.12 ± 2.25 


53309 


1 


29.83 ± 6.66 


26.00 ± 6.17 


7.00 ± 


3.78 


62.83 ± 8.99 


7.78 ± 


1.71 


6.78 ± 1.61 


1.83 ± 0.99 


16.38 ± 2.34 


16.38 ± 2.34 


53344 


4 


21.92 ± 6.77 


25.92 ± 6.17 


3.00 ± 


2.94 


60.83 ± 8.19 


5.81 ± 


1.63 


6.87 ± 1.64 


0.80 ± 0.78 


13.47 ± 2.17 


13.47 ± 2.17 


53345 


8 




















14.63 ± 0.57 


16.70-24.86 


53366 


3 




















8.78 ± 0.67 


9.69-15.45 


53366 


8 




















8.72 ± 0.67 


9.66-15.19 


53398 


1 




















11.31 ± 0.80 


12.44-19.78 


53398 


8 




















11.35 ± 0.66 


12.67-19.62 


53422 


7 




















8.23 ± 0.44 


9.23-14.16 


53741 


8 


14.00 ± 4.84 


12.00 ± 4.67 


3.92 ± 


3.18 


29.92 ± 6.55 


3.54 ± 


1.22 


3.04 ± 1.16 


0.99 ± 0.80 


7.67 ± 1.66 


7.57 ± 1.66 


53881 


2 


4.00 ± 3.18 


2.00 ± 2.66 


1.00 ± 


2.32 


7.00 ± 3.78 


3.26 ± 


2.69 


1.63 ± 2.17 


0.81 ± 1.89 


6.70 ± 3.08 


5.70 ± 3.08 


53891 


3 




















3.36 ± 0.43 


1.46-6.76 


53895 


7 


4.00 ± 3.18 


7.00 ± 3.78 


0.00 ± 


1.87 


11.00 ± 4.43 


2.06 ± 


1.64 


3.61 ± 1.95 


0.00 ± 0.96 


6.67 ± 2.28 


5.67 ± 2.28 


53947 





15.00 ± 4.97 


4.00 ± 3.18 


2.00 ± 


2.66 


21.00 ± 5.66 


3.77 ± 


1.26 


1.01 ± 0.80 


0.50 ± 0.67 


6.28 ± 1.42 


5.28 ± 1.42 


54002 


8 


24.67 ± 6.08 


14.00 ± 4.84 


2.92 ± 


2.94 


41.68 ± 7.58 


6.02 ± 


1.48 


3.42 ± 1.18 


0.71 ± 0.72 


10.15 ± 1.85 


10.15 ± 1.85 


54008 


9 




















10.64 ± 0.74 


6.52-13.24 


54052 


3 




















3.83 ± 0.47 


2.22-6.00 


54073 


9 


35.50 ± 7.07 


10.75 ± 4.43 


0.00 ± 


1.87 


46.08 ± 7.92 


1.55 ± 


0.31 


0.47 ± 0.19 


0.00 ± 0.08 


2.01 ± 0.34 


2.01 ± 0.34 


54114 


9 




















< 2.49 


< 2.49 


54145 


1 




















< 2.32 


< 2.32 


54169 


2 




















< 1.92 


< 1.92 
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Table 2 — Continued 





Soft 


Medium 




Hard 




Total 




SCR 






MCR 






HCR 




TCR 


ACIS-I Rate 


OBSMJD 


Counts 


Counts 




Counts 




Counts 


10~ 


3s 




1 


10~ 


3s 




1 


10~ 


3s 




1 


1A-3 -1 
10 s 


in-3 -1 
lU s 


54170.6 
















































< 0.88 


< 1.13 


54253.9 


6.00 ± 3.40 


2.92 ± 2.94 


1 


.00 


± 


2, 


.32 


8.92 ± 4.12 


1 


.02 


± 


0, 


,69 


0, 


,69 


± 





.60 


0. 


,20 


± 


0, 


,47 


1.82 ± 0.84 


1.82 ± 0.84 


54269.3 
















































< 3.28 


< 3.28 


54286.5 


4.00 ± 3.18 


1.00 ± 2.32 


0, 


.00 


± 


1, 


.87 


6.00 ± 3.40 


6, 


.12 


± 


4, 


,07 


1 


,28 


± 


2 


.97 


0. 


,00 


± 


2, 


.39 


6.40 ± 4.36 


6.40 ± 4.36 


54312.1 


9.00 ± 4.12 


6.00 ± 3.60 


1, 


.00 


± 


2, 


.32 


16.00 ± 6.09 


3, 


.47 


± 


1 


,59 


2, 


,32 


± 


1 


.39 


0. 


,39 


± 


0, 


.90 


6.18 ± 1.96 


6.18 ± 1.96 


54340.5 


6.83 ± 3.79 


9.00 ± 4.12 


0, 


,00 


± 


1 


.87 


15.83 ± 5.09 


1 


.70 


± 





,94 


2, 


,24 


± 


1 


.02 


0, 


,00 


± 





.46 


3.93 ± 1.26 


3.93 ± 1.26 


54369.6 


20.92 ± 5.66 


19.00 ± 5.44 


0, 


,00 


± 


1 


,87 


39.92 ± 7.38 


5, 


,29 


± 


1 


,43 


4, 


,81 


± 


1 


.38 


0, 


,00 


± 





.47 


10.10 ± 1.87 


10.10 ± 1.87 


54407.2 


18.00 ± 5.33 


9.00 ± 4.12 


1 


.83 


± 


2, 


.66 


28.83 ± 6.46 


4, 


,39 


± 


1 


,30 


2, 


,20 


± 


1 


.01 


0, 


,45 


± 





.65 


7.03 ± 1.58 


7.03 ± 1.58 


54411.6 
















































9.58 ± 0.71 


6.97-13.19 


54421.8 
















































6.30 ± 0.58 


4.45-8.81 


54431.2 


18.83 ± 5.45 


12.00 ± 4.57 


0, 


.83 


± 


2, 


.33 


31.67 ± 6.73 


4, 


,96 


± 


1, 


,44 


3, 


,16 


± 


1 


.20 


0, 


,22 


± 





.61 


8.34 ± 1.77 


8.34 ± 1.77 


54432.8 
















































7.69 ± 0.63 


6.84-11.57 


54441.6 
















































8.93 ± 0.70 


7.96-13.39 


54451.5 
















































12.39 ± 0.80 


11.23-18.35 


54463.7 


8.83 ± 4.12 


7.92 ± 3.96 


1 


.00 


± 


2, 


.32 


17.75 ± 5.33 


4, 


,23 


± 


1 


,97 


3, 


,79 


± 


1 


.90 


0, 


,48 


± 


1 


.11 


8.50 ± 2.55 


8.50 ± 2.55 


54617.5 


8.92 ± 4.12 


7.92 ± 3.96 


1 


.00 


± 


1 


.32 


17.83 ± 6.33 


4, 


,38 


± 


2, 


,02 


3, 


,89 


± 


1 


.94 


0, 


,49 


± 


0, 


.65 


8.75 ± 2.62 


8.76 ± 2.62 


54841.0 
















































< 1.96 


< 1.96 


54869.4 


















































a 


54901.6 


3.92 ± 3.18 


1.00 ± 2.32 


0, 


.00 


± 


1 


.87 


4.92 ± 3.40 


1 


,66 


± 


1, 


,26 


0, 


,40 


± 





.92 


0, 


,00 


± 


0, 


.74 


1.95 ± 1.35 


1.95 ± 1.35 


54980.1 
















































< 1.50 


< 1.50 



Note. — Columns 2, 3, and 4 arc the counts in the Soft (0.1-1.1 kcV), Medium (1.1-2 keV), and Hard (2-7 keV) band^^ respectively. The fourth 
column is Total Counts for ACIS observations (in the energy range of 0.1-7 keV). The next four columns are Soft, Medium, Hard, and Total count 
rates, respectively. All rates are corrected by a vignetting factor. The last column is the TCR converted into units of ACIS-I Count Rate, in the 
energy range of 0.1-7 keV (see text). 

^No upper limit could be calculated for rl-25 in this observation, as there were zero counts in both the source and background apertures. 



Table 3. Spectral Fitting Results 







Nh = 1 X 10^1 


cm 




Nh = 6.4 X 10^1 cm- 


-2 


Obs ID 


kT (keV) 


Normalization 


Lx (10^'^ ergs s'^) 


kT (keV) 


Normalization (10~®) 


Lx (10^^ ergs s"!) 


1575 
4720 
4721 
4722 


0.13 ± 0.01 
0.27 ± 0.04 
0.23 ± 0.02 
0.26 ± 0.02 


n ■7C+0.7 

^•^'-0.2 
9 1-1-0.3 
^•^-0.3 

2 0+°-3 
^■'-'-0.3 


0.4 
1.0 
1.1 
1.1 


0.08 ± 0.002 
0.19 d= 0.02 
0.16 ± 0.01 

U.10_o Q]^ 


s.o+\i 

14.01-1 
19.8l^i 


19.7 
4.2 
6.8 
9.1 



Note. — Spectral fits were carried out in the energy range 0.3-8.0 keV. Luminosities are given for the same energy range and 
are corrected for assumed absorption. All quoted uncertainties are la. 
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Table 4. Results of spectral simulations 



Sim obs ID 


Orig obs ID 


Nh (10^1 cm-2) 


SCR 


MCR 


HCR 


TCR 


4720 


1575 


1.1 


5.87 ± 2.42 


1.30 ± 1.14 


0.00 ± 0.04 


7.17 ± 2.72 


4720 


1575 


6.4 


7.62 ± 2.65 


0.99 ± 1.01 


0.00 ± 0.00 


8.60 ± 2.85 


4721 


1575 


1.1 


5.89 ± 2.44 


1.20 ± 1.09 


0.00 ± 0.05 


7.09 ± 2.70 


4721 


1575 


6.4 


7.51 ± 2.70 


0.95 ± 0.94 


0.00 ± 0.03 


8.46 ± 2.89 


4722 


1575 


1.1 


5.49 ± 2.47 


1.20 ± 1.10 


0.00 ± 0.06 


6.70 ± 2.67 


4722 


1575 


6.4 


7.24 ± 2.68 


0.92 ± 0.95 


0.00 ± 0.00 


8.16 ± 2.81 


1575 


4720 


1.1 


542.89 ± 23.59 


459.62 ± 21.86 


33.21 ± 5.65 


1035.72 ± 32.17 


1575 


4720 


6.4 


357.65 ± 19.76 


472.02 ± 21.16 


12.38 ± 3.66 


842.05 ± 30.02 


1575 


4721 


1.1 


720.86 ± 26.35 


430.37 ± 20.03 


17.61 ± 4.13 


1168.83 ± 33.24 


1575 


4721 


6.4 


531.22 ± 22.64 


484.71 ± 22.07 


6.41 ± 2.55 


1022.33 ± 31.11 


1575 


4722 


1.1 


655.97 ± 26.10 


484.68 ± 22.16 


27.87 ± 5.10 


1168.53 ± 35.45 


1575 


4722 


6.4 


618.62 ± 24.65 


477.95 ± 22.10 


4.58 ± 2.14 


1101.15 ± 32.47 



Note. — Col (1): Obs ID from which response files were taken from observation. Col. (2): Obs ID of model used for 
simulation (see Table 3). Col. (3): Nh value simulated. Cols (4), (5) and (6): Simulated total counts and associated Ict 
errors. 
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Table 5. pimms Method of Spectral Fitting for rl-25 in ACIS-S Observation 1854 



kT 


SCRp 


MCRp 


HCRp 




keV 


TCRp 


TCRp 


TCRp 


ergs s"-*^ 


0.025 


l.OOE+00 


4.55E-12 


O.OOE+00 


2.09E+39 


0.034 


l.OOE+00 


1.91E-08 


1.58E-19 


3.75E+38 


0.045 


l.OOE+00 


6.50E-06 


5.13E-14 


8.15E+37 


0.061 


l.OOE+00 


3.70E-04 


4.50E-10 


2.37E+37 


0.083 


9.94E-01 


5.99E-03 


2.94E-07 


9.28E+36 


0.112 


9.60E-01 


4.00E-02 


3.04E-05 


4.62E+36 


0.151 


8.64E-01 


1.35E-01 


7.55E-04 


2.84E+36 


0.203 


7.08E-01 


2.85E-01 


6.91E-03 


2.09E+36 


0.274 


5.32E-01 


4.35E-01 


3.29E-02 


1.81E+36 


0.370 


3.78E-01 


5.24E-01 


9.79E-02 


1.82E+36 


0.499 


2.61E-01 


5.30E-01 


2.09E-01 


2.09E+36 


0.673 


1.77E-01 


4.70E-01 


3.53E-01 


2.62E+36 


0.908 


1.19E-01 


3.83E-01 


4.99E-01 


3.52E+36 


1.224 


8.20E-02 


3.01E-01 


6.17E-01 


4.97E+36 



Note. — The spectral analysis in this tabic demonstrates the initial pimms method. We 
input various values oi kT (shown in Column 1) into pimms (assuming Nh = 1-1 x 10^^ 
cm^^) and calculated the ratios of the band count rates (shown in Columns 2-4). The last 
column shows the pimms prediction for the X-ray luminosity (see text). To determine which, 
if any, of these temperatures fit the source, we compare the pimms ratios to the observed ones. 
The observed ratios for this source are: 0.626 < fgf < 1.12, 2.65 x 10"^ < fg| < 0.262, 
-8.43 X 10-2 < < 8.43 x 10"^. The rows in bold are the models that fit the source in 
this observation. The range ol kT that fit was further reduced by breaking the soft band 
into smaller bins (see text). 
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Table 6. Binned Soft Band for rl-25 



ObsID 


MJD 


0.1-0.3 keV 


0.3-0.5 keV 


0.5-0.7 keV 


0.7-0.9 keV 


0.9-1.1 keV 


312 


51783.8 





2 


2 


3 





1581 


51891.1 











3 


3 


1854 


51922.4 


2 


8 


4 


5 


1 


1583 


52070.8 








1 


4 


2 


1575 


52187 


5 


34 


60 


46 


27 


4678 


52952.3 








2 


4 


2 


4679 


52969.9 








3 


1 


3 


4680 


53000.3 





1 





5 


3 


4682 


53148.7 








2 


6 


5 


4719 


53203.9 





2 


4 


7 


6 


4720 


53250.6 








7 


8 


10 


4721 


53282.9 





1 


10 


9 


9 


4722 


53309.1 





2 


8 


14 


6 


4723 


53344.4 








9 


6 


7 


7136 


53741.8 








1 


6 


7 


7137 


53881.2 








1 


2 


1 


7138 


53895.7 





1 


1 


1 


1 


7139 


53947.0 





2 


1 


6 


6 


7140 


54002.8 





2 


7 


10 


6 


7064 


54073.9 





2 


7 


13 


14 


( UUo 




n 

U 


n 

U 


1 
1 


Q 
O 


1 


8192 


54286.5 








2 


1 


1 


8193 


54312.1 





2 


1 


3 


3 


8194 


54340.5 








1 


4 


2 


8195 


54369.6 





1 


3 


8 


9 


8186 


54407.2 





2 


4 


5 


7 


8187 


54431.2 





2 


5 


5 


7 


9520 


54463.7 





1 


2 


3 


3 


9529 


54617.5 





2 


3 


3 


1 


10553 


54901.6 





1 





2 


1 
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Note. — Columns 3-7 are the raw counts in the 0.1-0.3, 0.3-0.5, 0.5-0.7, 0.7-0.9, and 
0.9-1.1 keV bands, respectively. 
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Table 7. Comparison of Results From XSPEC and pimms 







XSPEC 


pimms 


OBS ID 


OBSMJD 


kT (eV) 


kT (eV) 


1575 


52187.0 


120 - 140 


120 - 135 


4720 


53250.6 


230 - 310 


220 - 300 


4721 


53282.9 


210 - 250 


205 - 290 


4722 


53309.1 


240 - 280 


No Fit 



Note. — Results are for = 1.1 x 10^^ cm ^. The kT range for XSPEC are based on 
Icr uncertainties. For ObsID 4722, pimms was unable to determine a black body fit. 
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Table 8. rl-25 Spectrum 



ObsID OBSMJD kT (eV) LOG [Lx (ergs s-^)] 



312 


51783.8 


75 - 


- 110 


36.66 - 


37.24 


1581 


51891.1 


155 


- 200 


36.20 - 


36.32 


1854 


51922.4 


75 


- 83 


36.96 - 


37.09 


1583 


52070.8 


95 - 


- 140 


36.35 - 


36.89 


1575^ 


52187.0 


120 


- 140 


36.64 - 


36.69 


4678 


52952.3 


85 ■ 


- 185 


36.44 - 


37.40 


4679 


52969.9 


70 ■ 


- 135 


36.67 - 


37.71 


4680 


53000.3 


85 ■ 


- 185 


36.49 - 


37.48 


4682 


53148.7 


120 


- 205 


37.31 - 


37.83 


4719 


53203.9 


160 


- 205 


36.87 - 


37.06 


4720*^ 


53250.6 


230 


- 310 


36.90 - 


37.20 


4721* 


53282.9 


210 


- 250 


37.00 - 


37.20 


4722* 


53309.1 


240 


- 280 


37.00 - 


37.20 


4723 


53344.4 


200 


- 305 


36.90 - 


37.08 


7136 


53741.8 


245 


- 260 


36.71 - 


36.72 


7137 


53881.2 


95 ■ 


- 290 


36.54 - 


37.54 


7138 


53895.7 


175 


- 185 


36.76 - 


36.83 


7139 


53947.0 


95 - 


- 145 


36.95 - 


37.50 


7140 


54002.8 


160 


- 205 


36.94 - 


37.14 


7064 


54073.9 


120 


- 165 


36.40 - 


36.75 


7068 


54253.9 


100 


- 230 


36.17- 


36.96 


8192 


54286.5 


65 ■ 


- 145 


37.01 - 


38.25 


8193 


54312.1 


170 


- 315 


36.54 - 


36.84 


8194 


54340.5 


170 


- 250 


36.40 - 


36.64 


8195 


54369.6 


170 


- 265 


36.79 - 


37.05 


8186 


54407.2 


140 


- 160 


36.95 - 


37.09 


8187 


54431.2 


170 


- 230 


36.77- 


36.97 


9520 


54463.7 


180 


- 255 


36.73 - 


36.93 


9529 


54617.5 


190 


- 350 


36.68 - 


36.90 


10553 


54901.6 


65 - 


- 140 


36.55 - 


37.76 
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Note. — Results are for Nh = 1-1 x 10^^ cm^^. The minimum kT corresponds to the 
maximum luminosity for all rl-25 observations. Lx is in the 0.1-7 keV range. See text for 
how the kT and Lx were calculated. 

'^The spectra for these observations are from XSPEC fits (see text). The Lx for these 
observations were converted to the 0.1-7 keV range for consistency. 
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Table 9. rl-25 Color Ratios 



Detection Number 


OBSMJD 


Hard Color ((H-M)/T) 


Soft Color ((M-S)/T) 


1 


52187.0 


-0.13 


-0.74 


2 


53203.9 


-0.52 


0.03 


3 


53250.6 


-0.60 


0.34 


4 


53282.9 


-0.53 


0.20 


5 


53309.1 


-0.38 


0.10 


6 


53344.4 


-0.53 


0.24 


7 


53741.8 


-0.40 


0.20 


8 


54002.8 


-0.34 


-0.11 


9 


54073.9 


-0.38 


-0.24 


10 


54369.6 


-0.63 


0.25 


11 


54407.2 


-0.35 


-0.10 


12 


54431.2 


-0.51 


0.10 



Note. — The numbers in the first column correspond to the points plotted in Figure 5. 
Bands used for these ratios were: S:0.3-l keV, M:l-2 keV, H:2-8 keV, and T:0.3-8 keV. See 
text for how the Hard and Soft ratios were calculated. 



Table 10. ACS Observations of rl-25 



Date 


Dataset 


Filter 


Exposure Time (s) 


2004-01-23 07:14:05 


j8vp03010 


F435W 


2200 


2004-08-15 00:19:24 


j8vp05010 


F435W 


2200 


2006-02-10 03:58:08 


j9ju01010 


F435W 


4360 


2007-01-10 23:06:03 


j9ju06010 


F435W 


4672 



